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We use the functional renormalization group method to analyze the phase diagram of a 4-band 
model for the iron-pnictides subject to band interactions with certain Ai g momentum dependence. 
We determine the parameter regimes where an extended s-wave pairing instability with and without 
nodes emerges. On the electron-doped side, the parameter regime in which a nodal gap appears is 
found to be much narrower than recently predicted in On the hole-doped side, the extended s- 
wave pairing never becomes nodal: above a critical strength of the intra-band repulsion, the system 
favors an exotic extended d- wave instability on the enlarged hole pockets at much lower T c . At half 
filling, we find that a strong momentum dependence of inter-band pair hopping yields an extended 
s-wave instability instead of spin-density wave (SDW) ordering. These results demonstrate that an 
interaction anisotropy around the Fermi surfaces generally leads to a pronounced sensitivity of the 
pairing state on the system parameters. 

PACS numbers: 74.20.Mn, 74.20.Rp, 74.25.Jb, 74.72. Jb 



Introduction — The discovery of high-temperature su- 
perconductivity in iron arsenide and related compounds 
at the beginning of 2008 has triggered an enormous 
interest in condensed matter physics. This new class 
of materials exhibits transition temperatures T c beyond 
the conventional BCS regime upon electron 0] and hole 
doping 0] of a collinear antiferromagnetic parent state, 
with T c 's extending up to 56 K 0, thereby break- 
ing the cuprate monopoly on high-temperature super- 
conductivity 0- Present experimental evidence accom- 
panied by theoretical modelling suggest that the pair- 
ing in the iron-pnictides is different from the d-wave in 
cupratcs. This could be due to both the different strength 
of the electron-electron interactions in the two materials, 
which results in an itinerant antiferromagnet in the par- 
ent pnictide compounds, and to the topology of the Fermi 
surfaces (FS) with complex multi-band character [1, 
However, very similar to the cuprates is the belief that 
the magnetism of the parent state crucially influences the 
pairing symmetry of the doped system. 

Various approaches have been pursued to investigate 
the pairing symmetry in the iron-pnictides. By now, after 
a short period of analysis providing a wide-spread range 
of possible pairings [TO, 11 , [HJ > the general theoretical 
view started to converge to an extended s-wave order pa- 
rameter (denoted or s x 2 y 2) that takes opposite signs 
on the electron and hole pockets along the multi-band FS, 
which is consistent with some experimental data and also 
has broad theoretical support [i, HI, 14, T^, 1(| 17, 18| . 
However, experimentally, there is still no broad consen- 
sus about the nature of the order parameter. While most 
experiments can be explained in the framework of an 
gap 15[, several facts, such as the T 3 dependence of the 
NMR penetration depth over a significant temperature 
range, as well as the linear penetration depth in LaOFcP 



remain unsettled. Other experiments such as penetration 
depth on the 1111 and 122 compounds, as well as thermal 
conductivity can be explained by an order parame- 
ter but with large gap anisotropy. Gaps with significant 
momentum dependence, but no nodes, are reported ex- 
perimentally in [l9|. In contrast, ARPES data reveals 
very iso trop ic nodeless gaps on the hole Fermi surfaces 
j20l . |2ll [23, of magnitudes matching a strong-coupling 
form A(fc) = Aocos(fca;) • cos(k y ) [13[ in the unfolded 
Brillouin zone. Non-nodal gaps have been found in sev- 
eral analytical and numerical theoretical approaches to 
the pnictides [3, M, 0, EH 
ization group (fRG) studies 



26(| . Functional rcnormal- 
16| with orbital interactions 
reveal the presence of largely anisotropic, almost nodal, 
gaps. A recent RPA analysis suggests the existence of a 
nodal extended s-wave state within the 5-orbital Hub- 
bard model 



271 ]. It was also recently predicted that 



the momentum dependence of the interaction along the 
electron pockets may result in the development of gap 
anisotropy and possibly nodes on top of the constant s^ 1 



signa l I2S. I29H . In the strong coupling mean-field pic- 



ture [l3l. 15l| . the gap anisotropy is doping dependent: 
the gap has a form cos(fe x ) ■ cos(k y ) which becomes more 
anisotropic as the doping is increased. In a weak coupling 
expansion of FS interactions, the gap anisotropy can arise 
from the presence of an A\ g term cos(/ca;)-|-cos(A:y ) (which 
does not break the crystal symmetry but can create nodes 
on the (n, 0) and (0, 7r) electron surfaces) in the band in- 
teractions |28( | . 

In this Letter, we study the influence of interac- 
tion anisotropy on the superconducting instability by 
(fRG) [10, [H], methods. Our aim is to investigate the 
possibility that details of the material and the theoreti- 
cal model might play a decisive role in determining the 
pairing symmetry. In order to have a simple control and 
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FIG. 1: (Color online) (a) The band structure of the four 
band model. Electron pocket bands are denoted by /3i_2 and 
hole pocket bands by ai,2. The dashed line shows the Fermi 
level at half filling, (b) Change of Fermi surfaces upon doping. 
Bold lines are FS's at half filling. Upon electron (hole) doping 
(x — — .15 and x = .15 shown as dashed lines), the electron 
(hole) pockets grow while the hole (electron) pockets shrink, 
(c) Illustration of the patching segments in the Brillouin zone. 
The patches 1, ... ,24 and 49, . . . , 72 belong to the electron 
pockets situated at the M point, 25, . . . , 48 and 73, . . . , 96 to 
the hole pockets situated at the F point. 
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where i denotes the band index i = a±,a2, Pi, and 
the e"s are the on-site energies. The hopping dispersion 
parameters for the a bands around the T point are given 
by i£ = tf (cos k x + cos k y ) + t% cos k x cos k y , whereas pa- 
rameters for the (3 bands around the M point are given 
by = t± (cos k x + cos fc^) + if cos k x /2 cos k y /2. In units 
of eV and grouped by (e\ t\, i|), the parameters take on 
the values (-0.60,0.30,0.24) for a x , (-0.40,0.20,0.24) 
for a 2 , (1.70, 1.14, 0.74) for /3 U and (1.70, 1.14, -0.64) for 
/?2- The chemical potential can be adjusted such that two 
hole pockets emerge at the T point, whereas two electron 
pockets emerge at the M point. At half filling, i.e., /i = 0, 
hole and electron pockets are almost perfectly nested. 

Band interactions — We distinguish four types of band 
interactions. First, there are two types of inter-band in- 
teraction vertices depending on whether momentum from 
ingoing and outgoing propagators is transferred within 
the same band or between the bands: 



minimal model of the interaction anisotropy that allows 
to generate gap nodes, we do not consider interactions de- 
fined in an orbital representation, but use constant band 
interactions that only depend on whether the external 
legs are on hole- or electron-pockets. To these we add an 
Ai g momentum dependence in the inter-band pair hop- 
ping interactions qr with zero total momentum, as sug- 
gested recently in [281 ] . For the electron doped regime, 
we mostly find non- nodal s^, where the scale of the gap 
anisotropy increases with enhanced intra-band repulsion 
174. For highly dominant g<±, we indeed find a small pa- 
rameter regime of nodal s ± . On the hole-doped side, 
we first find a rather isotropic s signal where the scale 
of anisotropy does not decisively depend on 94. Upon 
further increasing (74 and hole doping, a phase transi- 
tion occurs, where both the electron and the hole pock- 
ets become nodal: the hole pockets develop a d-wave 
intra-band cooper pairing. This clearly demonstrates the 
possibility of a doping-dependent gap function. Further- 
more, we find that dominant momentum dependence of 
33 can increase T c of s^-wave pairing at half filling such 
that it even exceeds SDW as the previously leading in- 
stability. Recent experimental evidence for a supercon- 
ducting parent compound leaves open the possibility to 
have a leading superconducting instability even at half 
filling [31]. 

4-band model — As a tight-binding model, we constrain 
ourselves to the 4-band model given by Korshunov and 
Eremin [24[ containing the decisive weight of the den- 
sity of states in the vicinity of the FS's (sec Fig [lj) . In 
the folded BZ with two Fe ions per unit cell, the kinetic 
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with implicit sums over a, a' , and the band indices a and 
P extending over cvi^ and Pi,2, respectively (this conven- 
tion is kept throughout the article). While g x turns out 
to be rather unimportant with respect to the analysis 
of the leading instabilities driven by band interactions, 
g2 is necessary (but not sufficient) to drive the SDW in- 
stability in the (tt, tt) channel between electron and hole 
pockets in the folded BZ. Furthermore, we consider the 
inter-band pair hopping interaction 173: 

Hr= E 9i(d{ +qa A'- qa ,A^ d ^+^ (3) 

k.k'.q 

where, as one central point of our analysis, we include a 
momentum dependence of the Aig-projectcd pair hop- 
ping amplitude for the zero momentum cooper chan- 
nel [H: 

9 b s\k'=-k = 93 (1 + H cos k x + cos ky)) 

■ (1 + b(cos(k x + q x ) + cos(k y + q y ))) , (4) 

and constant 33 otherwise, b, the anisotropy scale, gives 
the relative scale of momentum dependence. To make 
connection to [l[, 53 defined in ^ is given in terms of 
the unfolded (u) momenta, which relate to the folded 
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(f) ones by k nx ,y = (fefii k{ y )/2. Finally, there is the 
intra-band pair interaction 
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where i extends over all band indices. For .ffj? an in the to- 
tal zero momentum Cooper channel, one Cooper pair be- 
longs to the electron pockets and the other one to the hole 
pockets, rendering g% to significantly deviate depending 
on the momenta along the electron FS, but gives only a 
constant value g% « 1 + 2b on the hole pockets. 

fRG method — We use the fRG method to study 
the flow of the two-particle coupling function 
V (fei, a, fc 2 , b, k 3 , c, &4, d), where ki, 2 (&3,4) denote 
the ingoing (outgoing) particles, a, . . . , d arc the dif- 
ferent band indices, and A is the energy cutoff above 
which the high energy contributions are integrated out 
and incorporated into the effective coupling function. 
Details on the implementation for the multi-band case 
of pnictides can be reviewed in While A; 4 is 

implicitly given by momentum conservation, the spin 
convention is chosen such that fei (k 2 ) and k 3 (k^) 
have the same spin. As usual, we omit the frequency 
dependence of the vertex, and assume that the relevant 
processes are located in the close vicinity of the FS. 
This, together with the neglect of self-energy corrections, 
is tantamount to a weak-coupling approach initially, 
although interactions strengths grow large under RG 
flow. To solve the RG equations numerically, the 
momenta in the BZ are discretized as shown in Fig. [1] 
Each momentum is confined within two neighboring 
dotted lines running from (jr, tt) to one of the corner 
points (n ± 7T, 7r ± 71"). For a given continuous momentum 
k and band index i, its value is effectively projected onto 
the Fermi momentum value of the ith band in the angle 
section where k is situated. Testing discretizations of 48, 
64, 96, 128, and 256 patches, we find that N = 96 both 
provides sufficient resolution and suitable computation 
time performance for our studies. Essentially, the RG 
flow in terms of the cutoff parameter A starts at cutoff 
energy scales of the order of the bandwidth, and succes- 
sively decreases towards the FS. In each differential step, 
both particle-particle and particle-hole contributions 
influence the evolution of the coupling function. Quite 
generically, these flows lead to strong coupling, i.e. one 
or more channels in the coupling function flow to large 
absolute values at a critical scale A c . Comparing the 
growth of the various pairing and density-wave channels 
allows one to compare between leading and subleading 
instabilities. For standard cases, the critical scale A c 
provides a reasonable estimate of the actual critical 
temperature T c for phase transitions when long-range 
ordering is possible. 

Electron doping — First we want to specify the pa- 
rameter window where nodal and non-nodal appears 
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FIG. 2: (Color online) (al)-(cl) Plot of the vertex function vs. 
incoming momenta for gi — 0. 1,(72 = 0.25,<?3 = 0.3,34 = 733, 
b = 1.1, and 7 varied as 2.4, 3, and 3.5 from left to right. k<t 
is fixed by total momentum, £3 is chosen to reside on patch 
sector 55. Upon increase of 7, one clearly observes the nodal 
structure on the electron pockets to become more pronounced, 
while the hole pockets show a homogeneously diverging vertex 
without anisotropy. (a2)-(c2) Plot of the associated supercon- 
ducting form factor. While the signal along the hole pockets is 
constant, the nodal peaks successively develop on the electron 
pockets. 



for electron doping (x = 0.2). We consider the case of 
</4 = 0.1 and (72 = 0.25, and choose (73 = 0.3 to be of 
the order of 32 (throughout the article, the interaction 
couplings arc given in units of eV where the total band- 
width is ~ 6eV). The 173 anisotropy scale b and 7 = 34/53 
span the relevant parameter space of interactions and al- 
low for nodal gaps according to Ref. [l|. For electron 
doping, the pockets nesting is lifted and the SDW insta- 
bility is suppressed. For 6^1, one finds a constant s^ 1 
instability, as observed in (23[, whose critical divergence 
scale T c decreases with increasing 7. The s ± instabil- 
ity manifests itself as a Cooper instability corresponding 
to divergent vertex couplings for fei = — k 2 , with a sign 
change of the vertex going from electron to hole pockets, 
as it can be seen both in the vertex plot and the form 
factors as shown in Fig. O Upon increasing b, the crit- 
ical divergence scale increases considerably and helps to 
counteract the intra-band repulsion g^. In addition, a gap 
variation starts to emerge on the electron pockets, while 
the hole pockets remain unchanged. In particular, the 
gap anisotropy increases upon increasing 7, which shows 
that the system favors a nodal variation to compensate 
the increasing intra-band repulsion 34 , as predicted in [l| . 
However, for 7^3, this anisotropy never becomes com- 
parable to the constant gap scale on the electron pockets, 
i.e. the nodes do not completely develop. Upon increas- 
ing 7 > 3, the gap variation on the electron pockets gets 
more pronounced, to finally yield a nodal extended s- 
wave instability (Fig.[2j;). However, for increasing b even 
further at constant 7, the value of 173 « 1 + 26 on the hole 
pockets is considerably enhanced and leads to a reduc- 
tion of gap anisotropy developing on the electron pock- 
ets. Hence, while the trends pointed out in Ref. [1( are 
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FIG. 3: (Color online) (al)-(cl) plot of the vertex function 
with inter-band interaction g\ = 0.1, <?2 = 0.25, gz = 0.3, 
54 = 753, 7 = 2.4, b = 1.1, x = -0.18, and 7 = 1.4, 2.0, 2.4. 
For (al) and (b2) we observe an ordinary a instability with 
small gap variation as shown in the form factor (a2) and (b2) . 
In (cl), we observe a leading signal on the hole pockets origi- 
nating from a d-wave pairing, where the nodes can be seen in 
the form factor (c2). (c3) gives a visualization of the extended 
d-wave state in the unfolded Brillouin zone. The bold dashed 
lines building the rhombohedron centered around the T points 
correspond to cos k x + cos k y yielding the sign change on the 
electron pockets around (±7r, 0) and (0, ±7r). The dashed lines 
along the diagonals of the BZ correspond to the d-wave func- 
tion cos k x — cos k y labelling the nodal points on the hole 
pockets. Upon folding, we observe that the gaps on the line 
r — » M have opposite signs. 

clearly visible, we only observe a clean nodal insta- 
bility with sign-changes around the electron pockets in a 
comparably narrow window of intermediate b and high 7 
(see Fig. [5] and Fig. O. In our understanding, this differ- 
ence occurs due to the renormalization of the interactions 
at higher energy scales before the pairing instability sets 
in. This effect is mentioned but not explicitly taken into 
account in and hence the parameter ranges for the 
nodal state differ. Concluding this part, we note that 
on the electron-doped side, the fully-gapped s±-state is 
rather stable. It is well possible that fine-tuning of the 
band structure like a strengthening of the scattering be- 
tween the electron pockets can enhance the anisotropy 
again, but primary instability favors an isotropic gap. 

Hole Doping — We now consider hole doping to the sys- 
tem, and choose the same interaction parameters g\ = 
0.1, 92 = 0.25, 53 = 0.3. For small anisotropy b, the 
system likewise develops a rather constant ,s ± instability, 
where T c decreases with 7. For considerable values of b 
and the ratio 7 as for electron doping, we find that the 
hole doped scenario still favors a constant instability 
(Fig. [3]). Unlike in the electron-doped case, increasing 
7 does not immediately induce gap anisotropy on the 
electron pockets: the only result (mostly) is a decreasing 
gap on the hole Fermi surfaces. This behavior can be ex- 
plained due to the small electron FS size: for hole-doping, 
the hole FS's dominate the behavior of the system. The 
Ai g term induces anisotropy mostly on the electron FS's, 
which at hole doping arc rather small. On the hole pock- 
ets, the A\ g term increases the size of the constant part 
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FIG. 4: (Color online) (al)-(cl) Plot of the vertex function 
with interaction couplings gi = 0.1, 32 = 0.25,p3 = 0.3,<?4 = 
0.85, a; = 0.0 and b varying from left to right as 1.5, 1.7, and 
1.8. Upon increasing b, the previously subleading s Cooper 
instability becomes more dominant and overcomes the SDW 
instability at b ~ 1.7. (a2)-c(2) s and SDW susceptibilities 
plotted for the different parameter settings. For b = 1.5, one 
observes the SDW susceptibility (dashed line) to be dominant 
w. r. t. to the s susceptibility. At b = 1.7, the divergence 
scales are nearly equal, and s becomes dominant for b = 1.8. 



of the interaction g^, which favors an s^. This, coupled 
to the fact that the electron FS's play a rather secondary 
role with respect to their hole counterparts, renders the 
symmetry to be for moderate 7. Increased 7 has the 
effect of reducing the scale of the superconducting insta- 
bility. Regarding the gap anisotropy, while in the elec- 
tron doped-case, increasing 7 resulted in nodal electron 
FS's (while keeping the hole FS's isotropic nodeless), for 
hole doping, the situation is completely different. Be- 
yond a critical 7, when the gaps have vanished, the 
system exhibits a phase transition with nodal supercon- 
ductivity on both hole and electron FS's. For 7 > 2.4, at 
comparably low transition temperature, we find that the 
leading instability becomes a d-wave intra-band Cooper 
pairing on the hole pockets. This is plausible from the 
FS topology upon doping as shown in Fig. [T]d: The hole 
FS's around T grow significantly upon hole doping. For 
dominant intra-band repulsion 34, the favorable ordering 
becomes a d-wave Cooper pairing on the hole pockets, 
as shown in Fig. [3] In terms of the hole FS only, this is 
a situation similar to the cuprate superconductors. We 
also find a subdominant d-wave type signal on the elec- 
tron pockets, which relates to the d-wave signal on the 
hole pockets shifted by tt so that there is still an overall 
sign change in the gaps going from the hole pocket to the 
electron pocket in the unfolded BZ; in this way, the sys- 
tem minimizes both intra and inter-band repulsion: this 
is the d-wave equivalent of the s instability, which one 
may denote extended d^-wave (Fig. (3^). This gap sym- 
metry represents another, time-reversal symmetric way 
to cope with the frustration pointed out in Rcf. 



34 1. 



While interesting in its own right, the phase diagram 
on the hole-doped side does not allow for a nodal electron 
FS while keeping the hole FS fully gapped. If this latter 
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FIG. 5: (Color online). Phase diagrams for different in- 
teraction settings, (a) 31 — 32 — 33 = 34 = 0.4 and (b) 
gi = 0.1, 32 = 0.25, 33 = 0.3, 34 = 733 = 0.65. The three axes 
are given by the filling factor x, momentum anisotropy b and 
the critical temperature T c . (a) In a dome around half filling, 
the SDW instability is leading, until it is overcome by the su- 
perconducting instability at higher b and electron (hole) dop- 
ing. T c increases with b, in particular on the electron doped 
side, (b) 7 is increased. The SDW dome around half filling 
shrinks in doping width but increases on the anisotropy line 
as the superconducting T c is decreased due to 34. For hole 
doping at comparably small amount of anisotropy b, we ob- 
serve nodal d . s gap variation on the electron doped side 
is enhanced, while true nodal extended s-wave is successively 
observed on the electron pockets only for even larger 7 ~ 3 
(see Fig. [6)1. 

gap structure turns out to be experimentally correct, we 
conclude that an A\ g symmetry term alone is not suf- 
ficient to make the electron FS nodal while keeping the 
hole FS fully gapped at both electron and hole doping. 

Interplay of SDW and at half filling — For all pa- 
rameter settings discussed in the previous sections, we 
always find a leading SDW instability as we approach 
half filling. This is because the SDW instability bene- 
fits enormously from the increased nesting of hole and 
electron pockets. However, we observe that upon fur- 
ther increasing b, the instability can overcome SDW 
and becomes the leading instability even at half filling. 
This is illustrated in Fig. [4] There, we find that upon 
changing b from 1.5 to 1.8, the s ± instability becomes 
dominant. The b c where this occurs increases with 7. 
Interestingly, the critical scale is not suppressed between 
the two regimes which in a less approximate treatment 
would bc most likely separated by a first order transition. 
This means that a slight change of the system parame- 
ters can turn the system from a high-scale SDW into a s ± 
state with comparable pairing scale - another possibility 
for a marked material-dependence of the phase diagrams 
of different pnictides, which may relate to recent mea- 
surements reported on in [33j |. 

Full Phase Diagram — For various interactions set- 
tings, we computed the complete phase diagram for dif- 
ferent fillings and anisotropy 6, four representatives of 
which are shown with increasing 7 in Fig. [5] and Fig. [6] 
We observe that for comparably small 7, the phase dia- 
gram looks rather uniform and only contains SDW and 
s with moderate gap variation as leading instabilities. 
For larger 7 and generally reduced T c due to strong intra- 
band repulsion, the phase diagram shows a more complex 




FIG. 6: (Color online) Phase diagrams for very strong intra- 
band repulsion 34 = 733. 31 = 0.1, 32 = 0.25,33 = 0.3,34 = 
733 = 0.80 (34 = 0.99) for (a) and (b), respectively, (a) Nar- 
row regime around half filling with SDW leading instability. 
T c is comparably small for the different phases, and, for small 
b and electron doping, undergoes our maximum resolution 
scale denoted by white parcels. The extended d wave regime 
on the hole doped side increases to larger b. (b) A small pa- 
rameter window with nodal s and small T c denoted by black 
is found for electron doping. 

structure including the extended d-wave instability on 
the hole doped side and the increasingly pronounced s 
gap variation for the electron doped side, which finally 
yields a small parameter window of nodal s for very 
large 7 (Fig. [BJa). Thus, we find smoothly connected in- 
stabilities at comparable scales, pointing to competing 
orders. This, however, suggests a potentially high sensi- 
tivity to material-specific parameters. 

Conclusion — We studied the effect of certain momen- 
tum dependence of band interactions for a 4-band model 
of the pnictides by adding an Ai g -symmetric term to the 
pair-hopping interactions, while keeping all other band 
interactions constant. While most of our data shows, in 
agreement with previous results, an wave supercon- 
ducting instability, on the electron doped side, we find 
that increasing intra-band repulsion enhances the gap 
anisotropy on the electron pockets, which may ultimately, 
but not readily, lead to a true nodal electron FS gap. For 
the hole-doped regime, however, the gap anisotropy re- 
mains rather small even on the electron FS upon increas- 
ing the interaction anisotropy, until a critical value of 
the intra-band repulsion beyond which the system favors 
an interesting state with intra-band d-wave cooper pair- 
ing on the hole pockets and a reminiscent extended and 
sign-reversed <i-wave signal on the electron pockets. This 
represents a novel way to reduce the repulsion within 
and between the Fermi surface pockets. However, in the 
range of our model, we arc unable to find a regime which 
satisfies the conditions of having gapped hole FS and 
nodal electron FS for both electron and hole doping. Fi- 
nally we demonstrated that even in the undoped state, 
small variations of the interaction parameters can turn 
the SDW ground state into a superconducting state with 
relatively high pairing scale. Taken together, the out- 
come of our studies is twofold. One possibility is that the 
bare interactions have little anisotropy around the Fermi 
surfaces and the main repulsive nesting is between elec- 
tron and hole pockets, then the phase diagrams should 
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feature fully-gapped s ± superconducting states. Alterna- 
tively, if the experiments confirm gap nodes, this would 
indicate a strong anisotropy of the bare interaction and a 
stronger role of other repulsive interaction, e.g., between 
the electron pockets. In this case, electron- and hole- 
doping may well lead to different pairing states, and even 
undoped superconductivity could be obtained for rela- 
tively mild parameter changes. One can ask how much 
these findings depend on the type of anisotropy of the 
bare interaction. In our understanding any pronounced 
wave vector-dependence of the interaction will very likely 
cause differences between electron and hole doping, as 
the relevant hole- or electron Fermi surfaces are located 
in different regions of the Brillouin zone. Likewise, for 
the undoped system, a strong interaction anisotropy gen- 
erally enhances pairing states competing with the SDW 
order such that material-specific parameter differences 
might lead to observable consequences even in the un- 
doped state. The trends observed in our study should 
hold for any specific choice of anisotropy in general. 
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